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Abstract 


Many numerical methods have been proposed in the last 30 years for inverse problems. 
While very successful in many cases, progress has lagged in other areas of applications which 
are forced to rely on limited-aperture measurements. In this paper, we introduce some tech- 
niques to retrieve the other data that can not be measured directly. We consider the inverse 
acoustic scattering of time harmonic plane waves and take the scattering amplitude to be the 
measurements. Assume that the scattering amplitude can only be measured with observation 
directions restricted in So, which is compactly supported in the unit sphere. Based on 
the reciprocity relation of the scattering amplitude, we prove a special symmetric structure 
of the corresponding multi-static response matrix. This will also be verified by numerical 
examples. Combining this, with the help of the Green’s formula for the scattered field, we 
introduce an iterative scheme to retrieve approximate full-aperture scattering amplitude. As 
an application, using a recently proposed direct sampling method [28], we consider the fast 
and robust sampling methods with limited-aperture measurements. Some numerical simula- 
tions are conducted with noisy data, and the results will further verify the effectiveness and 
robustness of the proposed data retrieval method and of the sampling method for inverse 
acoustic scattering problems. 


Keywords: Acoustic scattering, scattering amplitude, Multi-Static Response matrix, 
limited-aperture problem, direct sampling method. 
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1 Introduction 


The field of inverse scattering theory has been a large and fast-developing area of applied math- 
ematics for the past thirty years. The aim of research is to detect and identify the unknown 
objects through the use of acoustic, electromagnetic, or elastic waves. In the past thirty years, 
many numerical methods have been proposed to solve such kinds of problems , such as itera- 
tive methods, decomposition methods, linear sampling method, factorization method and direct 
sampling methods; we refer to [5, 6, 9, 10, 11, 13, 15, 17, 19, 20, 21, 22, 24, 26, 27, 28, 38, 39] 
and the references therein for these methods and some other related developments. The reader 
can also consult the recent monographs and review papers |3, 7, 8, 18, 32, 37] for a survey on 
the numerical methods. Most of the above research has considered full-aperture inverse scat- 
tering problems, i.e., the observation directions span the unit sphere. However, in many cases 


of practical interest, it is impossible to measure the data in all directions around the object- 
s, e.g., underground mineral prospection, mines locating in the battlefield, and anti-submarine 
detection. Correspondingly, in many studies|1, 4, 8, 14, 16, 18, 23, 25, 34, 35, 40], various re- 
construction algorithms have considered limited-aperture inverse scattering problems. Even if 
uniqueness of the inverse problems can be proved (see e.g., in [12, 29]), as one would expect, the 
quality of the reconstructions decreases drastically for this so called limited-aperture problem, and 
will actually deteriorates as the aperture decreases. Indeed, limited-aperture data can present 
a severe challenge for the numerical methods. A typical feature of the limited-aperture results 
is that the "shadow region" is highly elongated down range |8, 23]. Physically, the information 
available from the "shadow region" is very weak, in particular for high frequency waves|34]. For 
the two-dimensional problems the numerical experiments of Decomposition Methods in [16, 40] 
indicate that satisfactory reconstructions need an aperture not smaller than 180 degrees. 

All the above research use the limited-aperture data directly for the corresponding inverse 
problems. Our main contribution in this paper is, based on the model, to introduce some 
techniques to retrieve the data that can not be measured directly. Thus the limited-aperture 
problem has actually been changed into the classical full-aperture problem. We take as our model 
problem the acoustic scattering by time-harmonic plane waves. In this work, the measurements 
are only taken over a limited range of angles. Contrary to this, as the first step, we use plane 
waves from all directions. 

We begin with the formulations of the acoustic scattering problems. Let k = w/c > 0 be the 
wave number of a time harmonic wave where w > 0 and c > 0 denote the frequency and sound 
speed, respectively. Let Q C R"(n = 2,3) be an open and bounded domain with Lipschitz- 
boundary 00 such that the exterior R"\Q is connected. Furthermore, let the incident field u be 
a plane wave of the form 


u(x) = u(x; d) = et, xz eER”, (1.1) 


where d € S$" denotes the direction of the incident wave and S"~1 := {x € R” : |x| = 1} is the 
unit sphere in R”. Then the scattering problem for the inhomogeneous medium is to find the 
total field u = u’ + uê such that 


Aut+tk(1+qu=0 in R”, (1.2) 
n—1 ðu! 
li "2 | — — iku |) = 1. 
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where q € L®(R”) such that S(q) > 0 and q = 0 in R”\Q, the Sommerfeld radiating condition 
(1.3) holds uniformly with respect to all directions ĉ := x/|z| € S"~+. If the scatterer Q is 
impenetrable, the direct scattering is to find the total field u = ut + uê such that 


Aut k?u=0 in R"\Q, (1.4) 
B(u) =0 on dQ, (1.5) 
dim or E -— itu’) = 0, (1.6) 


where B denotes one of the following three boundary conditions: 


(1) B(u) := u on OQ; (2) B(u) := au on OQ; (3) Ba) = - + àu on 02 


corresponding, respectively, to the case when the scatterer 2 is sound-soft, sound-hard, and of 
impedance type. Here, v is the unit outward normal to 0Q and A € L®(ƏN) is the (complex 
valued) impedance function such that S(A) > 0 almost everywhere on OQ. Uniqueness of the 
scattering problems (1.2)—(1.3) and (1.4)-(1.6) can be shown with the help of Green’s theorem, 
Rellich’s lemma and unique continuation principle, see e.g., [12]. The proof of existence can be 
done by variational approaches (cf. [12, 36] for the Dirichlet boundary condition and |7, 33] for 
other boundary conditions) or by integral equation methods (cf.[12, 30, 31)). 

Every radiating solution of the Helmholtz equation has the following asymptotic behavior at 
infinity |18, 28]: 


. n—2 . 
iF e ikr 1 
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uniformly with respect to all directions ĉ := x/|z| € 9”71. The complex valued function u®(#) = 
u% (ĉ; d) defined on the unit sphere $"~! is known as the scattering amplitude or far-field pattern 
with ĉ € $”~! denoting the observation direction. 

It is well known that the scatterer Q can be uniquely determined by the scattering amplitude 
u®(#,d) for all ĉ,d € S"~! [12]. Due to analyticity, for uniqueness it suffices to known the scat- 
tering amplitude on a subset oo C $”—! with nonempty interior. Unfortunately, the analytic 
continuation is a strongly ill-posed problem, which is established by the diabolical theorem (see 
Atkinson [2]). To the author’s knowledge, there exists no sufficient numerical method for the 
analytic continuation in our case. 

This paper is organized as follows. In section 2, we begin with an introduction of the multi- 
static response (MSR) matrix, which is the scattering amplitude in the finite case. With the 
help of the reciprocity relation of the scattering amplitude, we find that the MSR matrix can 
be regarded as a 2-by-2 block matrix with special symmetric properties. Based on this fact, we 
are able to retrieve directly and exactly part of the scattering amplitude in special angles. In 
subsection 2.2, using the Green’s formula, we introduce a technique to compute the scattering 
amplitude in the rest angles. Together with the technique proposed in subsection 2.1, a novel 
algorithm is proposed to retrieve full-aperture data from limited-aperture data. As an application, 
in section 3, we introduce some direct sampling indicator functionals for inverse acoustic problems 
by using only limited-aperture data or retrieved full-aperture data. Some numerical simulations 
in two dimensions will be presented in section 4 to verify our novel algorithms. 

For simplicity, in later sections, we restrict our presentation to the two-dimensional case. 
The three-dimensional analysis, although more complex, presents essentially no further compli- 
cations. Actually, all the results of this paper remain valid in three dimensions after appropriate 
modifications of the fundamental solution and the radiation condition. 


2 Date retrieval from the model 
In RÊ, we choose an equidistant set of knots 6; := (i —1)m/m, i = 1,2,...,2m from [0, 27). 
Assume that we have a set of incident plane waves with incident directions 
di := (cos 0;, sin 0;), i = 1,2,..., 2m. 
The scattering amplitude are measured in different observation directions 


ĉj := (cos 0j, sin 04), j = 1,2,..., 2m. 


In the finite case we define the multi-static response (MSR) matrix F puu € C?Mx2?M by 


CO CO CO 
uig Ute ‘` Ulam 
CO CO CO 
U31 U22 "°° U2 2m 
F fuu = : ” ; (2.1) 
Co O CO 
Ugm,1 U2m,2 °° U2m,2m 


where uyg = u™ (ĉj; di) for 1 < i,j < 2m corresponding to 2m observation directions ĉj and 
2m incident directions d;. The MSR matrix F fuu given in (2.1) is regarded as the scattering 
amplitude in full-aperture. 

In practical applications, the scattering amplitude can only be measured in a limited-aperture. 
After a necessary rotation of the coordinate axes, we may make the assumption without loss of 


generality that we can take the first / columns of F fuy to obtain 


0O oO [0,6] 
ui Ute 07° Uja 
OO [0,6] OO 
D 2,1 2,2 2 
Finit = : , 1<1<2m. (2.2) 
(e e] OO OO 
U2m,1 U2m,2 °° ° Yam 


(1) E€ C2mxt 


In other words, a realistic inverse problem is to reconstruct Q from the matrix Finit 


2.1 Structure of the MSR matrix and its application 


We want to remark that F fuu is NOT a symmetric matrix, i.e., F fuu £ Full’ Here and through- 
out the paper we use the superscript ”T” to denote the transposition of a matrix. We can 
partition both the rows and columns of the 2m-by-2m MSR matrix Ffu to obtain a 2-by-2 
block matrix 


n) 
F ful = ( a ) (2.3) 


21 22 


where F;; € C™*™,i,7 = 1,2 designates the (i,j) block (submatrix). The following theorem 
indicates the structures for each submatrix. 


Theorem 2.1. Fi; = F}, Fig = Fi, and Fo = F}. 


Proof. Recall that the scattering amplitude is unchanged if the direction of the incident field and 
the observation direction are interchanged |12], i.e., 


u™®(ĉ, d) = u®(—d,—ĉ), for all ĉ,d € St. (2.4) 


For all urs € Fin, using the reciprocity relation (2.4), we have 


, 


WG = wd) 
= u™(-di;—ĉ;) 
= u™(-—(cos6;, sin 6;); —(cos 0j, sin 0;)) 
= u™((cos(b; + 7), sin(6; + 7)); (cos(@; + 7), sin(0; + 7))) 
= u™((cos Oi4m, Sin him); (cos Oj4m, Sin Aj4m)) 
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= u% 1<i,j<m. 


j+m,i +m) 


Thus, we have Fi; = F4. 
Similarly, For all Urjim © F12, using the reciprocity relation (2.4) again, we have 


oo = 
Uij+m = 


cos(6; + 7), sin(O; + 7)); (cos(0j;4m + 7), sin(0;4m + 7))) 
cos im, Sin him); (cos 8;, sin 6; )) 


= tag. [Sheen 


( 
a 


Thus, we have F12 = F3. The equality Fo; = F} can be treated analogously. 


Clearly, the above results also hold for phaseless MSR, matrix. If we set 


A R F 
12 Fu 

F ful = | i : ' 
22 Fə 


then F putt is symmetric, i.e., F put = Fip We want to emphasize again that the MSR matrix 


F fuu is NOT symmetric. 

As a direct consequence of Theorem 2.1 one has immediately the following technique. 

First technique to retrieve new data. All the data a 1<i<2m,l <j < 2m are 
missed in the measurements. However, by using Theorem 2.1, the following data 


O O 
Li “m+1, u 


co . 
m 2 m+1,2m 
o0 oo oo 
=(l) en m+2,I+1 “m+2,1+2 m+2,2m 
FD a = E i i , 1<1<2m, (2.5) 
oo oo oo 
Umi “m+, +2 °° Um-+1,2m 


have actually also be exactly retrieved. Here, we have set URS E U amj if i > 2m,1 < j < 2m. 


2.2 Green’s formula and its application 


Let B be a bounded domain with connected complement such that Q C B and the boundary 
OB is of class C?. Let v denote the unit normal vector to the boundary OB directed into the 
exterior of B. We recall that the fundamental solution ®(z, y), x € R?,x Æ y, of the Helmholtz 
equation is given by 


®(0,y) = THY) (kle — y), (2.6) 


where HP is the Hankel function of the first kind of order zero. Disregarding the scattering 


objects, the scattered field u'(-;d) is a radiating solution to the Helmholtz equation in R?\B. 
Then we have Green’s formula |12] 


T = usl: PEEL) — PE YD 54 g E 
wed =f fema Gow Daey) dsg), ceRAB (27) 


Letting x tend to the boundary OB, with the help of jump relations, it can be shown that (¢, Y) := 
(us a) loz € H'/?(@B) x H~‘/?(@B) must be a solution of the following two boundary integral 


Ov 
equations 
se) =2 f fom FEY — pajalan) base), © € a8, (28) 
B ð O® (x, y) 7 . p 7 
va) = 2z f {ow FOX -uee baw), ce ae. (2.9) 


For later use, we define the space 
W :={(¢,~) € H"? (3B) x H~-/?(OB) : ($, Y) is a solution to (2.8) — (2.9).}. 


Based on the Green’s formula (2.7), it is well known that the scattering amplitude u(-; d) 
has the following form (cf. [18]) 


—ikê- s 

u” (ĉ; d) = a {wud — iy aeren) ds(y), ĉeest. (2.10) 
An important observation is that the Cauchy data (u a) loz is independent of the variable 
&. From (2.10) we find that the scattering amplitude can be computed in any direction if the 
Cauchy data (e ae) loz is known in advance. Let St be the measurement surface, which is 
an open part of the unit sphere St with nonempty interior (open relative to $1). If we already 
know the scattering amplitude in S4, then it is natural to look for the Cauchy data (us, ou" | loz 


by solving the following integral equation 
F(9(-54), (5 4))(@) =u (â; d), FE Sp, (2.11) 


where F : W - L?(S}) is defined by 


en ike. D 
Fecawaya@ = | {ow Vode} dtu), tesh (12) 


ðB 
Theorem 2.2. The operator F : W — L? (SL) is compact, injective with dense range in L? (S$). 


Proof. The operator F is certainly compact since its kernel is analytic in both variables. 


Let ($, Y) € W satisfy F(¢(; d), Y(;d))(ĉ) = 0 in $j. By analytic we have 
Be tkt-y E 
oly; d) ——— — yly;d eter asty =0, ĉes! 2.13 
[foun - va w) (2.13) 
Note that the left hand side of (2.13) is actually the scattering amplitude of the scattered field 


wê given by 


s — P OP(x,y) _ : x s z 2D 
w= f foma REN —vuraeace.y} dst), ce RAB, 


By Rellich’s lemma, from (2.13) we obtain w° vanishes in R?\B. Now, using jump relations 
yields 


= £ eaa £ s P 
o=300)+ f {ow REN —vma(en}asy), zean, 


ae £ ? OE) £ s £ 
o= SH tos f {ow eX — van baw, ae aB. 


Recall that (¢, Y) € W, which implies (¢, y) is also a solution of the boundary integral equations 
(2.8)-(2.9). Hence, ¢ =~ = 0 and F is injective. 

We consider the adjoint F* of F and show that it is injective as well which proves the denseness 
of the range of F. For all h € L(,9$) we extend h by zero in $1\S4 to obtain h € L?(S'). Recall 
the Herglotz wave function vp of the form 


un(y) := L eY ®h(z)ds(@), y € R?. 


Then we obtain that the adjoint operator F* : L?($}) > H~'/?(0B) x H'/?(B) is given by 


ie fa: Ovp, 
F*h= (St -vs) ba (2.14) 


Indeed, by interchanging the order of integration, we have 


nike: pa 
Povie = f f {0m0 - wide) dyads) 


A fois DERE — yy: ama) as(y) 


= (6). (F-m) 


= ((6,~),F°h), 


where the last two equalities are understood in the sense of dual paring (H'/?(@B) x 
H-‘/?(0B), H~‘/?(0B) x H'/?(8B)). Here and in the following, Z denotes the complex con- 
jugate of z € C. 

We proceed by showing that the adjoint operator F* is injective. Let h € L? (SL) be such 
that F*h = 0 on OB. Again extending h by zero in S1\S4 to obtain h € L?(S1). We find that 
the Cauchy data of the Herglotz wave function v, vanishes on OB. Note that the Herglotz wave 
function vp is an entire solution of the Helmhlotz equation in R?. Thus by Holmgren’s uniqueness 
theorem we deduce that vp, vanishes identically in R. This further implies that h = 0 on OB 
[12] and the proof is finished. 


Second technique to retrieve new data. As described at the begin of this section, 
assume that we have 2m equidistantly distributed incident directions in S1, and we choose the 
first l directions as the observation directions in St. Then, practically, we obtain the scattering 


amplitude as shown in the matrix FO Theorem 2.2 indicated that, for every d; € S', we 
may find a pair of approximate solution (¢,~) of the equations (2.8)-(2.9) and (2.11) by using 
the Tikhonov regularization. Inserting this into formula (2.10), we then obtain the approximate 


scattering amplitude in other directions in S ae 


2.3 From limited-aperture to full-aperture 


Based on the two techniques proposed in the previous subsections, we introduce an algorithm 
to retrieve the scattering amplitudes that cannot be measured. 


Data retrieval scheme: 


(1) 


e Step 1. Collect the limited-aperture scattering amplitude F 7), , that be measured directly. 


Set s = 0, which is the iterative numbers. 


(2) 


e Step 2. Based on the first technique, we obtain directly the data FO a shown in (2.2). 


e Step 3. Based on the second technique, compute the scattering amplitude in the unavailable 
observation directions 


Mnrew := {ĉip Shia Êlpt, tom st t+1; 2m st—t+2,°°° a spacer he 


That is, we compute in the 2t new directions close to the known directions. Ifl +t < m, 
set l =l +t,s = s + 1 and return to Step 2. 


e Step 4. Using the first technique again to obtain the full-aperture scattering amplitude 
Fit- 


We believe the iterative method used in the scheme will increase the accuracy of the full- 
aperture scattering amplitude F)_,-. Clearly, the first technique gives more accurate data than 
the second one. However, we can only retrieve part of the data based on the first technique. Thus, 
a good approximate solution of the equations (2.8)-(2.9) and (2.11) from the second technique 
is highly desired. 

We finally remark that we make no use of the boundary conditions or topological properties 
of the underlying object Q. In other words, the full-aperture data is retrieved without knowing 
in advance the type of the boundary condition or the number of components. 


3 Applications to sampling methods for inverse scattering prob- 
lems 


Recently, in [28], we proposed a novel direct sampling method for inverse acoustic scattering 
problems based on the following indicator functional 


I(z) := |$(z; —d)F pun” (z; 2)”, (3.1) 


where b(2; —d) _ (ekz, ged oe ,e7 tkz-dam ) and olz; ĉ) = (ekzi eikzê2, a , ekz 22m), 
Clearly, only matrix multiplications are involved in the computation, thus it is very fast and 
robust against measurement noise from the numerical point of view. More importantly, the 
indicator is independent of any a priori information of the unknown scatterers. Theoretically, 
the indicator functional J (z) has a positive lower bound if the sampling point z located inside the 
scatterer, and decays like the bessel functions as the sampling points away from the boundary of 
the scatterer. Thus it is expect that the indicator takes its maximum on or near the boundary 
of the scatterer. 


A natural modification of the above sampling method for limited-aperture data is to introduce 
the following indicator 


l T a\12 
Diimie 2) = |(25 -DF ni bimi 2 ê), (3.2) 
where Quimit(z; 2) := (e8721 etkzê2 ... , e’***1) corresponds to the limited-aperture observation 


directions and F® is the limited-aperture data given by (2.2). 


limit 
Recall the limited-aperture data FO given in (2.5), which is nearly exactly reconstructed 


by the known data FO Define limit (z; —d) i= (ete dmt1 emikzdm+2 ... | e—the-dmti) and 
btimit (232) := (e87 e~ the Fit2 |... | e~thet2m) Then, based on the first data retrieval tech- 


nique introduced in subsection 2.1, one may also consider the second indicator 
l A T mf) Fi A 
Limit (2) = (2; -AJF lini Ptimiel( 25) + dtimie(25— A) Flim 0 (5 @)?. (3.3) 


We expect that the quality of the reconstructions can be improved by using the indicator J/,__,,.(z). 
Finally, based on the scheme introduced in the previous section, we also introduce the fol- 
lowing indicator 


Thun (z) == |9(2; —d)Fi+ ro" (z; ê). (3.4) 


As will seen in the next section, the quality of the reconstructions indeed improved greatly by 
using the indicator If,,)(2). 


4 Numerical examples and discussions 


Now we turn to present some numerical examples in two dimensions to illustrate the applicability 
and effectiveness of our data retrieval techniques and novel sampling methods for inverse prob- 
lems. All the programs in our experiments are written in Matlab and run on a Core i5-5200U 
2.2GHz PC. 

There are totally three groups of numerical tests to be considered, and they are respectively 
referred to as BlockSymmetric, DataRetrieval and SamplingMethods. 

For the numerical examples, we consider the scattering of plane waves by a cylinder with a 
non-convex kite-shaped cross section with boundary OQ illustrated in Figure 1 and described by 
the parametric representation 


Kite: z(t) = (cost + 0.65 cos 2t — 0.65,1.5sint), O<t< 2r. (4.1) 


A commonly used criterion for judging the quality of a reconstruction method is whether the 
concave part of the obstacle can be successfully recovered. The typical feature in the limited- 
aperture problems is that the concave part usually failed to be reconstructed if the the observation 
angles do not cover the concave part of the obstacle. For the synthetic experiments in this 
section, we used the boundary integral equation method to compute the scattering amplitudes 
Ung Pd = 1,2,:-- ,2m, for 2m equidistantly distributed incident directions and 2m observation 
directions. These data are then stored in the matrices F fun € C2mx2m_ We further perturb the 
scattering amplitude up”, by random noise using 


00,6 _ , 00 . 
Ung = Upg + O(r1 + ira), 
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True domain 


ol 


4b 


6 4 5 0 2 4 6 
Figure 1: Kite shaped domain for numerical examples. 


where rı and rg are two pseudo-random values drawn from the standard uniform distribution 
on the open interval (—1,1). The value of 6 used in our code thus presents the absolute error 
level, which is the amount of physical error in a measurement. Denote by Feu the full-aperture 


measurements. Let Re, be the matrix from the first l column of Fe ana and thus denote 
the limited-aperture data. We have used 4m equidistantly distributed nodes to discrete the 


boundary OQ of the scatterer. 


Example BlockSymmetric. This example is designed to verify Theorem 2.1. We set the 
wave number k = 1 and m = 8, i.e., the scattering amplitude is collected in 16 observation 
directions and 16 incident directions. The four block matrices are given as follows. It is easy to 
check that 


Fu =F, Fuo=Fip and Fo, =F}. 
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—1.6275 + 0.60221 —0.7591 + 0.77612 0.2507+0.55132 0.9295 + 0.37092 

ie, = | —1.1468 + 0.62822 —1.5928 + 0.51822 —0.8128+ 0.96107 0.4977 + 0.99792 
—0.4468 + 0.83447 —1.1669 + 0.3184i —1.5897 + 0.45867 —0.8128+0.9610i |” 
—0.0929 + 0.81012 —0.4366 + 0.58097 —1.1669 + 0.31842 —1.5928 + 0.51822 
—1.6275 + 0.60222 —1.1468 + 0.6282: —0.4468 + 0.8344; —0.0929 + 0.81012 

Fa = | —0.7591 + 0.77612 —1.5928+4 0.5182: —1.1669 + 0.3184: —0.4366 + 0.58092 
0.2507 + 0.55137 —0.8128 + 0.96107 —1.5897 + 0.45861 —1.1669 + 0.31847 |’ 
0.9295 + 0.37097 0.4977 + 0.99792 —0.8128 + 0.96102 —1.5928 + 0.51822 
1.3970 + 0.09502 0.9295 + 0.37092 0.2507 + 0.55137 —0.7591 + 0.77614 

Fo = | 0.9295 + 0.37092 0.4875 + 0.02132 0.5576 + 0.07802 0.1625 + 0.77582 
0.2507 + 0.55137 0.5576 + 0.07802 0.7922 + 0.32322 0.4381 + 0.97782 : 
—0.7591 + 0.77612 0.1625 + 0.77582 0.4381 + 0.9778: 0.1905 + 1.03721 
0.0471 + 0.71347  —0.0929+ 0.81012 —0.4468 + 0.83442 —1.1468 + 0.62822 

Fo = —0.0929 + 0.81017 0.19054 1.03722 0.4381+0.9778¢ 0.1625 + 0.77581 


—0.4468 + 0.83441 0.43814 0.97782 0.7922 + 0.382322 0.5576 + 0.07802 
—1.1468 + 0.62827 0.16254 0.77582 0.5576 + 0.07802 0.4875 + 0.02132 


Example DataRetrieval. In this example, we tested the validity of algorithm proposed in 
section 2. The same as the previous example, the underlying sound soft obstacle is a kite shaped 
domain as shown in (4.1). We took wavenumber k = 5, and compute the scattering amplitude 
u™(ĉ;, di), i,j = 1,2,--- ,256, for 256 incident and observation directions. Assume that the 
scattering amplitude can only be measured in the first 64 directions, i.e., 2; = (cos 0;, sin 6;) 
with 6; € [0, 7/2), j =1,2,--- ,64. The artificial domain B is chosen to be a ball centered at the 
origin with radius 5. We have used 64 equidistantly distributed nodes to discrete the boundary 
OB of the artificial domain B. 

Figure 2 shows the scattering amplitudes in all observation directions corresponding to two 
incident directions d = (1,0) and d = (0,1) with absolute error level 6 = 0.02. We note that 
the data carry a considerable amount of noise since F(1,37) = —0.0101 + 0.46497. As expected, 
the closer to the measurable directions, the better of the data reconstructions. Considering 
the severe ill-posedness of data retrieval for analytic functions |2], the reconstructed data can be 
regarded to be a good approximation of the full-aperture measurements. We have also considered 
the simplest and most common type of interpolation uses polynomials. Figure 2(a)-(b) show a 
polynomial p(x) of degree 16 that fits very well the first 64 data in a least squares sense. However, 
as shown in Figure 2(a)-(b), the method of polynomial interpolation indeed does not work at 
all for our data retrieval problems. The corresponding result with 6 = 0.1 is shown in Figure 3. 
Surprisingly, the imaginary part seems very robust to the noise. Such a fact is also observed in 
the later examples, see Figures 4(b,d) and 5(b,d). 

Figure 4 shows the result with incident direction d = (—0.9808, 0.1951). We observe that the 
scattering amplitude are reconstructed very well for the observation angles in [7, 27). Actually, by 
symmetric property of the block matrix F12, these data are the same as the scattering amplitude 
u™ (#49, di), i = 1,2,--- , 128. Here, the observation direction ĉ249 = (0.9808, —0.1951) is very 
close to measurement angles in [0,7/2), thus could be reconstructed well by using the second 
technique. 
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For incident angles in [r, 37/2], we have obtained nearly exact scattering amplitude for 
all observation directions. In Figure 5 we show results for two directions d = (—1,0) and 
d = (0,—1). This further verify the special symmetric structure of the multi-static response 
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Figure 2: Example DataRetrieval. Comparison between true data and reconstructed data 
with different incident directions d = (1,0), (0,1). The absolute error level 6 = 0.02. 


Example SamplingMethods. As an application, we consider the numerical methods for 
inverse acoustic scattering problems by using the limited-aperture data. In the simulations, we 
used a grid G of 121 x 121 equally spaced sampling points on some rectangle |[—6,6] x [—6, 6]. 
For each point z € G, we compute the indicator functionals given in (3.1)-(3.4). 

The resulting reconstruction by using the indicator functional [jjnj4(z) with limited-aperture 
en is shown in Figure 6(a). We see that the illuminated part is well 
constructed, but the shadow region is highly elongated down range. This is typical of limited- 


aperture results. As shown in Figure 6(b), the reconstruction is improved a little by using 


scattering amplitude F 
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Figure 3: Example DataRetrieval. 
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with different incident directions d = (1,0), (0,1). The absolute error level 6 = 0.1. 
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Figure 4: Example DataRetrieval. Comparison between true data and reconstructed data 
with different incident directions d = (—0.9808, 0.1951). (a)-(b), 6 = 0.02; (c)-(d), 6 = 0.1. 
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Figure 5: Example DataRetrieval. 
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with different incident directions d = (—1, 0), (0,—1). The absolute error level 6 = 0.1. 
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the indicator Jj,,,,,. Combining the retrieved scattering amplitude, the quality of the target 
reconstruction has been improved greatly by using Jj,,,,;,, and Ifun, respectively, as shown in 
Figure 6(b) and Figure 6(c). In particularly, by using the indicator functional If,(z), the two 
wings of the kite appear and the shadow region is reconstructed very well. As a comparison, the 
result by using the indicator functional I(z) with full-aperture scattering amplitude measurement 
Fy is shown in Figure 6(d). Considering the severe ill-posedness of the data reconstruction 
of an analytic function and the absolute error level 6 = 0.1, the target reconstruction given in 


Figure 6(c) is satisfactory. 


5 Concluding remarks 


The limited-aperture problems arise in various areas of practical applications such as radar, re- 
mote sensing, geophysics, and nondestructive testing. A typical feature of the numerical method 
with limited-aperture data is that the illuminated part can be reconstructed well, while the shad- 
ow domain failed to be recovered. In this paper, based on the mathematical model, we introduce 
some data retrieval techniques to approximately reconstruct the missing data that can not be 
measured directly. Both theoretical foundation and numerical experiments are presented. Using 
the retrieved full-aperture data, the direct sampling method proposed in a recent paper [28] yields 
satisfactory reconstruction of underlying objects. 
We conclude with some remarks for future work. 


e The second technique proposed in subsection 2.2 is to solve an ill-posed problem. We have 
used the Tikhonov regularization with regularization parameter a = 1071. As observed 
in Figures 2 and 3, the retrieved data get worse if the direction is far away from the 
measurable directions. A fast and stable method for solving the equations (2.8)-(2.9) and 
(2.11) is highly desired. 


e Recall the following integral equation of the first kind 
Soh = ur, (5.1) 


where the integral operator So : L?(0B) + L?(S") is defined by 


(S$..0)(@) = fe ¥6(y)dstu) (5.2) 
OB 

The integral equation (5.1) plays an important role in the Kirsch-Kress Decomposition 
method |19], where the auxiliary surface OB is chosen to be contained in the unknown 
scatterer Q. To retrieve full-aperture data, one may firstly solve the equation (5.1) by the 
Tikhonov regularization in L?($}) and then insert the solution ¢ into u%°(#) := (S.0@) (2) 
to obtain the missed data. Unfortunately, the integral operator S. has an analytic kernel 
and therefore equation (5.1) is severely ill-posed. Besides, one has to assume that B is 
contained in the unknown scatterer Q, which is not reasonable in many applications. We 
hope to be able to address this method and report the progress elsewhere in the future. 


e Of great practical importance would be limited-aperture data with not only observation 
directions, but also the incident directions. In particular, a realistic case is to consider the 
back-scattering limited-aperture data. 
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(b) Tiimit (z) 


(c) Ifuu(z) (d) I(z) 


Figure 6: Example SamplingMethods: Shape and location reconstructions by using different 
indicators. The absolute error level ô = 0.1. 
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e It would be interesting and useful to consider the scattering by buried objects, where the 
measurements are only available in the upper half space. The data for the inverse problem is 
available over a limited aperture, which implies that the solution of the inverse problem will 
be degraded compared to situations in which data can be gathered on a sphere containing 
the object. 
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